High-speed droplet migration in silicon 



H. E. Clifie and T. R. Anlhony 

General Elcaric Research and Development Center, Schenectady. New York 12301 
(Received 27 August 1975; in Trout form 6 February 1976) 

llijh-icmpermurc ihcrtromigtahon of aluminum-rich droplets in silicon has resulted in droplcl migration 
rales up to 10 cm/day. These unexpectedly high rates are an order of magnitude greater than those 
previously reported in the literature. The large enhancement of droplet migration rales at high temperatures 
is shown to be caused by a substantial increase in the flow of liquid across the droplet and by a favorable 
variation in the change of the composition of the liquid at high temperatures. The variation of droplet 
velocity with droplet size is reported over a large range of droplet sizes and temperatures for both the 
<100> and <lll> migration directions in silicon. The undcrsaturation and supersaturalion required to 
make a (1 1 1) silicon plane dissolve and deposit, respectively, are calculated from the data. Finally, from the 
velocity of large droplets, the liquid interdiffusion coefficient for the aluminum-silicon system is 
determined. 

[■ACS numbers: o6..10.Kv. 8l.20.lly 



I. INTRODUCTION 

Liquid droplets in a solid migrate in a thermal gradi- 
ent because atoms of the solid dissolve into the liquid 
at the hot interface of the droplet, diffuse across the 
droplet, and deposit on the cold interface of the drop- 
let. The resulting flux of dissolved solid atoms from the 
hot to the cold side of the droplet causes the droplet to 
migrate towards the hot end of the crystal. 

Previous experiments in a variety of materials 1-15 
have shown droplet migration rates of less than 1 cm/ 
day over a limited range of temperature. Since the 
transport of dissolved solid atoms across the droplet is 
limited ultimately by the magnitude of the liquid dif- 
fusion coefficient, high temperatures were not expected 




to markedly increase droplet migration rates because 
of the small temperature dependence of liquid diffusion 
coefficients. As a result, droplet migration experiments 
were never extended to relatively high temperatures. 

However, as is shown in this paper, the equation for 
droplet migration must be modified at high tempera- 
tures to Include both a variation in the change of com- 
position of the droplet liquid with temperature and a 
p!iy6icHi "iow of 'liquid that occurs across a droplet dur- 
ing migration. The resulting migration equation indi- 
cates a sharp increase in droplet migration rate at high 
temperatures. This prediction is confirmed experi- 
mentally in this paper by high-temperature droplet mi- 
gration rates that are an order of magnitude greater 

than those previously reported in the literature. 15 -" 

II. EXPERIMENTAL PROCEDURE 

Two phosphorus-doped lO-flcm, silicon single crys- 
tals containing 10 4 dislocations/cm 1 and with respective 
(100) ruid (111) crystallographic orientations along their 
cylindrical axis wore obtained from Texas Instruments 
Corp. From each 2.5-cm-diam by 40-cm-long crystal, 




. FIC. 1 . The mask used to etch a scries of square holes of ly deposited on the cooler face of the silicon wafer migrate up 

varying size In a silicon wafer by photolithography techniques. the temperature gradient through the wafer. 
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FIG. 7. The variation of droplet vclorld, wiik ^ „ i i ■ 
crease £ff ct ^.raUon^^^^^ . 

sprtir K e ?° rded WHh Pt " 10% Rh thermocouples in- 
serted , n holes in the side of the sample. After dronlet 
migration, the waiers were sllced J Mel 
of dr 0 p, e t s and the drop.et migration distance was de- 
termined from infrared transmission microphotographs 

CI OSS S«CttOn« Of th« rWrt »r ? .i»S left r-ohlnH » h „ ^ „„„ ' 

lets were then revealed by sectioning and stainimf fn 
several cases, the shapes of the solidified droplet were 
exaimned in the scanning electron microscope after 
the surrounding silicon had been etched away with I 
solution of 10 HNO,, 4 HAc, and 1 HF 



Sn^X?^",™" <1 ! 1) l '" rcctl - »- rovca.o.1 by the do- 
<«■ transmission photomicrograph 20x 

III. EXPERIMENTAL RESULTS 
A. Migration in the <1 11 ) direction 

luff. 3). Droplets smaller than 0.0175 cm did not ml 
pSem C, ' yStfll b ° CnUSe ° f a Surracc ^«Sr 
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aropiei niiKi-aii, 1(! i„u,<. g00) dlroctlon 
<l transmission photomicrograph along 
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the (100) direction as pyramids bounded by four for- 
ward (111) planes and a rear (100) plane (Fig. 8). The 
twisted doped trails loft behind some droplets migrating 
in thc<100> direction Indicate that the four (111) facets 
do not always dissolve at a uniform rate. Nonuniform 
dissolution of the four forward (111) facets can also 
cause the regular pyramid shape of the droplet to be- 
come distorted into a trapezoid as shown experimental- 
ly in Fig. 9 and schematically in Fig, 10. 

A scanning electron micrograph of a small droplet 
from which the surrounding silicon was selectively 
etched away is shown in Fig. 11. In this droplet, a 
4-/J-dcep ledge was seen on a dissolving (111) facet. 
At higher magnifications, kinks are observed on the 
ledge (Fig. 12). Evidence of such large dissolution 
ledges is also seen in infrared transmission of another 
droplet in Fig. 8. 

Like the (111) direction, the droplet migration rate 
along the (100) direction also decreases with decreasing 
droplet size (Fig. 13). A comparison of Fig. 13 with 
Fig. 6 indicates that the decrease in droplet velocity 
with decreasing size in the (100) direction is twice that 
in the (111) direction. 



IV. DISCUSSION 

A. Effect of the temperature on the migration rate 

To explain the marked increase in droplet velocity at 
high temperatures, a migration rate equation must be 
derived (o include the effects of the mass flow of fluid 
across a droplet during migration and the variation in 
the change of composition of the liquid with tempera- 
ture. Since the latter factor can be included in a 
straightforward way, let us concentrate for a moment 
on the effects of fluid flow across the droplet. 

We consider the case of a binary droplet migrating in 
a thermal gradient in a solid (aluminum-rich liquid 
droplet in silicon) where one component (aluminum in 
our experiment) has a negligibly small solid solubility 
in the solid. A coordinate system moving with the drop- 
let at the droplet velocity V is chosen for convenience 
since the droplet interfaces are stationary in this frame 
of reference. Conservation of mass of component 1 
(silicon) at the rear interface of the droplet requires 
that the flux J l of component 1 passing across the drop- 
let be equal to the removal rate VCf of component 1 at 
the rear interface nf the droplet. 



= VCf . 



(1) 



Here Cf is the concentration of component 1 in the solid 
phase. A similar mass conservation equation J^- VCf 
hnlris for component 2 (aluminum). However, since we 
are assuming that the solubility of component 2 is neg- 
ligible (Al in solid Si) in the solid, the mass conserva- 
tion equation for component 2 reduces to 



(2) 



in the frame of reference of the moving droplet. That is, 
alt of component 2 is carried along with the migrating 
droplet. 

At first glance, the lack of a flux of component 2 is a 



little be' --Rome since we know that a concentration 
gradient of con. r •-><mt I mid thus of component 2 Is im- 
posed across the dropiet by the temperature gradient. 
Consequently, we would expect a diffusion flux of com- 
ponent 2 similar to the diffusion flux of component 1 
which generates droplet migration. Although such a 
diffusion flux of component 2 does occur In the forward 
direction In the droplet, it Is exactly canceled out by a 
flow or drift of the liquid 16 - 11 towards the rear Interface 
of the droplet. If D 2 Is the intrinsic diffusion coeffi- 
cient of component 2 in the liquid and U' is the drirt 
velocity 16 -" of the tiquld In the droplet, with respect to 
the moving interface, we can express this cancellation 
between the diffusion flux and fluid flow for component 



Z as 



2 2 a* 2 



(3) 



For component 1 on the other hand, the fluxes gen- 
erated by the concentration gradient and the fluid flow 
are both towards the rear of the droplet and hence re- 
inforce each other, causing the droplet to migrate more 
rapidly. This fluid flow U< which Increases with in- 
creasing temperature [because of the decrease of C[ 
with increasing temperature In Eq. (3)1 is the direct 
cause of thr> marked inrrea"? in droplet rr.;gr::!:oa rate 
at high temperatures. The total flux of component 1 in 
the droplet produced both by pure diffusion and fluid 
flow is 



- D. — 



(4) 



(5) 



*-U'C'. 

The two components of the liquid satisfy the relation 

V t Cl+V t Ci =1, 

where V, and V, are respectively the partial molar vol- 
umes of components 1 and 2 in the liquid. 

Combining Eq. (1)— (5) to eliminate J u J z , t>€i/!>z, 
and W, the velocity V of the migrating droplet is found 

v= Vlc[lii + v[C[Dt 3C[ (6) 

By using the definition of the interdiffusion coefficient 
D = V' 2 C : l D l + V\C[D 2 , this expression simplifies to 

V —JL-*£L (7 , 

At low temperatures in the Al-Si system where the 
droplet is largely liquid aluminum Cjt' 2 '~l and Eq. (7) 



which has been used previously to interpret liquid drop- 
let migration in semiconductors. 5 " 6 ' 12 Here D t is the 
intrinsic diffusion coefficient of Si in liquid aluminum, 
Cf is the number of moles/unit volume of silicon \,\ 
solid silicon and dCjdz is the concentration gradient 
of the silicon in the liquid aluminum droplet. 

At higli temperatu/es on the other hand, the phase 
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diagram of the Al -Si system causes the liquid droplet 
to be composed almost entirely of Si and V[C[ approach- 
es unity, the velocity then becomes 



where D 2 is the intrinsic diffusion coefficient of Al in 
liquid Si. Since Cjl'j-^l, the droplet velocity becomes 
very large al high temperatures as we have observed. 
In the limit as C^-0, 'he velocity of the droplet does 
not become infinite however, since our assumptions of 
uniform fluxes and neglect of the latent hent evolved and 
absorbed respectively on the depositing and dissolving 
interfaces breakdown al extremely large droplet 
velocities. 

B. Interface kinetics and the interdiffusion coefficient 

The composition gradient SC|/3z in the liquid droplet 
can he calculated by dividing the concentration drop ACJ 
across the droplet by the droplet thickness £(dC{/9« 
= ACj/t). The concentration drop ACJ is generated 
by the charge in solubility (iC E /HT)CL induced by the 
temperature change CL across the droplet minus the 
undersatu ration U and supersatu ration S required to 



. Those factors 



droplet migration r: 



VELOCITY I0' 6 cm/uc 
FIG. 1 5. The sum or the unricrsnturntlon V nml supers.nu ra- 
tion S required to cause respectively silicon dissolution nnd 
deposition on the forward and renr faces of a droplet migrating 
In the (10n> direction at a rate appropriate to Urn droplet 
vccloclty V. 

cause dissolution and deposition on the forward and rear 
interfaces of the droplet lACJ = {SC E /dT)GL - U + S\.<- > 



3Cl 3C £ . 

~ = TT C 



u+s 



(8) 

In the Al-Si system, the temperature gradient C in the 
droplet is equal to the applied thermal gradient because 
the thermal conductivities of the solid and liquid are 
nearly equal."' 7 

7. The interdiffusion coefficient 

At large droplet sizes we may neglect the interface 
Wnelic term in Eq. 8. The droplet velocity from Eqs. 
(7) and (8) (hen becomes 



(9) 



which may be used to calculate the interdiffusion coef- 
ficient D as shown in Fig. 14. With D = D a exp( - Q/RT), 
values of Z3 0 = 5xiO-' cmVsec and <? = 3 kcal/mole were 
obtained which are in close agreement with liquid dif- 
fusion coefficients measured by other techniques."' 23 
Since the interdiffusion coefficient depends only weakly 
on temperature, it is apparent that the large variation 
in the droplet migration rate with temperature (Fig. 5) 
is caused by the large changes of \/C\v[ and 3C £ /3T 
with temperature (see Table I). 

2. Interface kinetics 

From the variation of droplet velocity V with droplet 
thickness L [see Eqs. (7) and (8)], the undersatu ration 
V and supersaturations S required respectively on thj 
dissolving r.r.rf depositing faces of the droplet can be 
calculated. For droplets migrating in the (100) 
direction, 
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U+S = 50K » t0 '» fractions 
cm/sec 

where V l00 is the velocity of the droplets In the <100) 
direction In cm/sec and U + S Is expressed In concen- 
tration units of atom fractions. Since both the pyramid 
droplets migrating along the (100) and the triangular 
platelets migrating along the (111) have (111) facets for 
their forward and rear faces, one would expect that the 
U + S values, alter a correction for the effective thick- 
ness L of the pyramid droplet and the resolved normal 
velocity of its (111) facets, would be equal. In fact, both 
the (100) data and the (111) data indicate that U + S (111 
facets) = 200V 1U (atom fraction) (cm/sec)' 1 , where V m 
is the normal resolved velocity of a 111 facet plane in 
cm/sec. 

C. Migration velocity 

A more useful expression for the migration velocity 
may be derived in terms of atomic fractir- >' of the 
droplet element (in our case aluminum). Y.;e equilibrium 
composition of the liquid is given on a phase diagram of 
(he system of interest. We convert our expression (9) 
to atomic fraction using the relations C l V l =X and 
0,1', = ! _,\\ We approximate {\/C' s ){dC E /hT) by the 
liquidus slope BX/ST because volume changes are 
usually small. We obtain an exuression for the velncitv 

>°M£> <>«> 

that clearly shows the effect of enchanced velocity as 
X~ 1 which occurs at high temperatures. The diffu- 
sivity of most liquids are similar and do not vary much 
with temperature. Also the value of the thermal gradient 
in the solid is not expected to vary much for a given 
apparatus. Thus the main factor determining the migra- 
tion velocity is the solubility in the liquid as a function 
of temperature. 

Doping semiconductors to produce electronic de- 
vices' 4-32 may be done much more rapidly by droplet 
migration than by conventional solid-state diffusion. A 
significant amount of time may be saved by using high 
temperatures to obtain high migration velocities. 
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